*******************************************************
*******************************************************
*******************************************************
*** CREATED BY JONATHAN KING AND JESSICA SCHOENHERR ***
***** REPLICATION DATA FOR "A MATTER OF OPINION?" *****
*******************************************************
*******************************************************
*******************************************************

use "AbortionData20220224.dta"

**********************
*** BASELINE MODEL ***
**********************

reg decisionTherm ib5.treatmentGroup
est store abortionBaseline

est restore abortionBaseline
margins, at(treatmentGroup=(1(1)5))

est restore abortionBaseline
margins, at(treatmentGroup=(5)) ///
	at(treatmentGroup=(1)) post
test _b[1._at] = _b[2._at]

est restore abortionBaseline
margins, at(treatmentGroup=(5)) ///
	at(treatmentGroup=(2)) post
test _b[1._at] = _b[2._at]

est restore abortionBaseline
margins, at(treatmentGroup=(5)) ///
	at(treatmentGroup=(3)) post
test _b[1._at] = _b[2._at]

est restore abortionBaseline
margins, at(treatmentGroup=(5)) ///
	at(treatmentGroup=(4)) post
test _b[1._at] = _b[2._at]

est restore abortionBaseline
margins, at(treatmentGroup=(1)) ///
	at(treatmentGroup=(2)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionBaseline
margins, at(treatmentGroup=(1)) ///
	at(treatmentGroup=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionBaseline
margins, at(treatmentGroup=(1)) ///
	at(treatmentGroup=(4)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionBaseline
margins, at(treatmentGroup=(2)) ///
	at(treatmentGroup=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionBaseline
margins, at(treatmentGroup=(2)) ///
	at(treatmentGroup=(4)) post
test _b[1._at] = _b[2._at]

est restore abortionBaseline
margins, at(treatmentGroup=(3)) ///
	at(treatmentGroup=(4)) post
test _b[1._at] = _b[2._at]
* not

********************
*** PARTISANSHIP ***
********************

* dem 1
* ind 2
* rep 3

reg decisionTherm ib5.treatmentGroup##ib1.partisanship
est store abortionPartisanship

est restore abortionPartisanship
margins, at(treatmentGroup=(1(1)5) partisanship=(1(1)3))

* democrats across treatment
* compared to treatment group 1
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(1)) ///
	at(treatmentGroup=(1) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(1)) ///
	at(treatmentGroup=(1) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(1)) ///
	at(treatmentGroup=(1) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(2) partisanship=(1)) ///
	at(treatmentGroup=(1) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to treatment group 2
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(1)) ///
	at(treatmentGroup=(2) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(1)) ///
	at(treatmentGroup=(2) partisanship=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(1)) ///
	at(treatmentGroup=(2) partisanship=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(1)) ///
	at(treatmentGroup=(3) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(1)) ///
	at(treatmentGroup=(3) partisanship=(1)) post
test _b[1._at] = _b[2._at]
*significant p=0.0369; 72.56 (3) vs 62.54 (4)

*compared to treatment group 4
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(1)) ///
	at(treatmentGroup=(4) partisanship=(1)) post
test _b[1._at] = _b[2._at]
* not 0.05 level, yes at 0.1 level p = 0.0512


* republicans across treatment
*compared to treatment group 1
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(3)) ///
	at(treatmentGroup=(1) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 58.76 (5) vs 44.73 (1), p = 0.0167

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(3)) ///
	at(treatmentGroup=(1) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(3)) ///
	at(treatmentGroup=(1) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(2) partisanship=(3)) ///
	at(treatmentGroup=(1) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(3)) ///
	at(treatmentGroup=(2) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(3)) ///
	at(treatmentGroup=(2) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(3)) ///
	at(treatmentGroup=(2) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(3)) ///
	at(treatmentGroup=(3) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(3)) ///
	at(treatmentGroup=(3) partisanship=(3)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(3)) ///
	at(treatmentGroup=(4) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* not

* independents across treatment
*compared to treatment group 1
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(2)) ///
	at(treatmentGroup=(1) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(2)) ///
	at(treatmentGroup=(1) partisanship=(2)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(2)) ///
	at(treatmentGroup=(1) partisanship=(2)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionPartisanship
margins, at(treatmentGroup=(2) partisanship=(2)) ///
	at(treatmentGroup=(1) partisanship=(2)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(2)) ///
	at(treatmentGroup=(2) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(2)) ///
	at(treatmentGroup=(2) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(2)) ///
	at(treatmentGroup=(2) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

*compared to treatment group 3
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(2)) ///
	at(treatmentGroup=(3) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(2)) ///
	at(treatmentGroup=(3) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

*compared to treatment group 4
est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(2)) ///
	at(treatmentGroup=(4) partisanship=(2)) post
test _b[1._at] = _b[2._at]
* not

* dems vs. reps within a treatment
est restore abortionPartisanship
margins, at(treatmentGroup=(1) partisanship=(1)) ///
	at(treatmentGroup=(1) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 67.5 (1) vs 44.7 (3), p = 0.000

est restore abortionPartisanship
margins, at(treatmentGroup=(2) partisanship=(1)) ///
	at(treatmentGroup=(2) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 64.7 (1) vs 52.1 (3), p = 0.0186

est restore abortionPartisanship
margins, at(treatmentGroup=(3) partisanship=(1)) ///
	at(treatmentGroup=(3) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 72.6 (1) vs 49.7 (3), p = 0.0000

est restore abortionPartisanship
margins, at(treatmentGroup=(4) partisanship=(1)) ///
	at(treatmentGroup=(4) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 62.5 (1) vs. 51.1 (3), p = 0.0273

est restore abortionPartisanship
margins, at(treatmentGroup=(5) partisanship=(1)) ///
	at(treatmentGroup=(5) partisanship=(3)) post
test _b[1._at] = _b[2._at]
* yes, 71.2 (1) vs 58.8 (3), p = 0.0222

**************
*** GENDER ***
**************

reg decisionTherm ib5.treatmentGroup##i.female
est store abortionGender

est restore abortionGender
margins, at(treatmentGroup=(1(1)5) female=(0(1)1))

* men across treatment
*compared to treatment group 1
est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, 66.2 (5) vs 52.6 (1), p = 0.0114

est restore abortionGender
margins, at(treatmentGroup=(4) female=(0)) ///
	at(treatmentGroup=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(3) female=(0)) ///
	at(treatmentGroup=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(2) female=(0)) ///
	at(treatmentGroup=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(2) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 66.2 (5) vs 56.3 (2), p = 0.0724

est restore abortionGender
margins, at(treatmentGroup=(4) female=(0)) ///
	at(treatmentGroup=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(3) female=(0)) ///
	at(treatmentGroup=(2) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 66.2 vs 58.0, p = 0.1229

est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(4) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, 66.2 (5) vs 55.1 (4), p = 0.0340

* women across treatment
*compared to treatment group 1
est restore abortionGender
margins, at(treatmentGroup=(5) female=(1)) ///
	at(treatmentGroup=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 65.6 (5) vs. 59.1 (1), p = 0.1530

est restore abortionGender
margins, at(treatmentGroup=(4) female=(1)) ///
	at(treatmentGroup=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(3) female=(1)) ///
	at(treatmentGroup=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(2) female=(1)) ///
	at(treatmentGroup=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionGender
margins, at(treatmentGroup=(5) female=(1)) ///
	at(treatmentGroup=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 65.6 (5) vs 58.2 (2), p = 0.1026

est restore abortionGender
margins, at(treatmentGroup=(4) female=(1)) ///
	at(treatmentGroup=(2) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionGender
margins, at(treatmentGroup=(3) female=(1)) ///
	at(treatmentGroup=(2) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionGender
margins, at(treatmentGroup=(5) female=(1)) ///
	at(treatmentGroup=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 65.6 (5) vs 61.9 (3), p = 0.4135

est restore abortionGender
margins, at(treatmentGroup=(4) female=(1)) ///
	at(treatmentGroup=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionGender
margins, at(treatmentGroup=(5) female=(1)) ///
	at(treatmentGroup=(4) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 65.6 (5) vs. 58.6 (4), p = 0.1302

* men vs women, within a treatment
est restore abortionGender
margins, at(treatmentGroup=(1) female=(0)) ///
	at(treatmentGroup=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 52.6 (0) vs 59.1 (1), p = 0.1823

est restore abortionGender
margins, at(treatmentGroup=(2) female=(0)) ///
	at(treatmentGroup=(2) female=(1)) post
test _b[1._at] = _b[2._at]
* no 56.3 (0) vs 58.2 (1), p = 0.6985

est restore abortionGender
margins, at(treatmentGroup=(3) female=(0)) ///
	at(treatmentGroup=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 58.0 (0) vs 61.9 (1), p = 0.4262

est restore abortionGender
margins, at(treatmentGroup=(4) female=(0)) ///
	at(treatmentGroup=(4) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 55.1 (0) vs 58.6 (1), p = 0.4683

est restore abortionGender
margins, at(treatmentGroup=(5) female=(0)) ///
	at(treatmentGroup=(5) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 66.2 vs 65.6, p = 0.9061

******************
*** FULL MODEL ***
******************

reg decisionTherm ib5.treatmentGroup##ib1.partisanship##i.female
est store abortionFull

est restore abortionFull
margins, at(treatmentGroup=(1(1)5) partisanship=(1(1)3) female=(0(1)1))

* democrat men across treatment
*compared to treatment group 1
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, 77.4 (5) vs 61.84 (1), p = 0.0415

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, 77.4 (5) vs 61.6 (2), p = 0.0427

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 77.4 (5) vs 68.4 (3), p = 0.2197

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(4) partisanship=(1) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, 77.4 (5) vs 57.1 (4), p = 0.0039

* democrat women across treatment
*compared to treatment group 1
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 67.4 (5) vs 70.6 (1), p = 0.5754

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* not

*compared to treatment group 2
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 67.4 (5) vs. 66.8 (2), p = 0.9208

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(3) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 67.4 (5) vs. 76.3 (3), p = 0.1638

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(3) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(4) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 67.4 (5) vs. 66.6 (4), p = 0.8838

* Republican men across treatment
*compared to treatment group 1
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 49.8 (5) vs. 44.8 (1), p = 0.5751

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* not

*compared to treatment group 2
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 49.8 (5) vs. 52.8 (2), p = 0.7372

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 49.8 (5) vs. 48.4 (3), p = 0.8749

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(4) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no, 49.8 (5) vs 52.5 (4), p = 0.7639

* Republican women across treatment
*compared to treatment group 1
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 64.5 (5) vs. 44.65 (1), p = 0.0127

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not 

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 2
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 64.5 (5) vs. 51.5 (2), p = 0.0950

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 3
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 64.5 (5) vs. 51.0 (3), p = 0.0722

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
*not

*compared to treatment group 4
est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(1)) ///
	at(treatmentGroup=(4) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 64.5 (5) vs. 49.5 (4), p = 0.0698

* democrat men vs women, within a treatment
est restore abortionFull
margins, at(treatmentGroup=(1) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(1) partisanship=(1))

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(2) partisanship=(1))

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(3) partisanship=(1))

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(4) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(4) partisanship=(1))

est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(5) partisanship=(1) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(5) partisanship=(1))

* republican men vs women, within a treatment
est restore abortionFull
margins, at(treatmentGroup=(1) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(1) partisanship=(3))

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(2) partisanship=(3))

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(3) partisanship=(3))

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(4) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(4) partisanship=(3))

est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(3) female=(0)) ///
	at(treatmentGroup=(5) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* not

est restore abortionFull
margins, dydx(female) at(treatmentGroup=(5) partisanship=(3))

* men, republican vs. democrat, within a treatment
est restore abortionFull
margins, at(treatmentGroup=(1) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0281

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(1) female=(0))

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(2) female=(0))

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0068

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(3) female=(0))

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(4) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* no

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(4) female=(0))

est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(0)) ///
	at(treatmentGroup=(5) partisanship=(3) female=(0)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0016

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(5) female=(0))

* women, republican vs. democrat, within a treatment
est restore abortionFull
margins, at(treatmentGroup=(1) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(1) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0002

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(1) female=(1))

est restore abortionFull
margins, at(treatmentGroup=(2) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(2) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0326

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(2) female=(1))

est restore abortionFull
margins, at(treatmentGroup=(3) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(3) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, p = 0.0003

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(3) female=(1))

est restore abortionFull
margins, at(treatmentGroup=(4) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(4) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* yes, 0.0200

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(4) female=(1))

est restore abortionFull
margins, at(treatmentGroup=(5) partisanship=(1) female=(1)) ///
	at(treatmentGroup=(5) partisanship=(3) female=(1)) post
test _b[1._at] = _b[2._at]
* no, 0.6786

est restore abortionFull
margins, dydx(partisanship) at(treatmentGroup=(5) female=(1))

***************************
*** WITH ADDED CONTROLS ***
***************************

* destring income, education, and age
replace income = "." if income == "NA"
destring(income), replace

replace education = "." if education == "NA"
destring(education), replace

reg decisionTherm ib5.treatmentGroup##ib1.partisanship##i.female income education age
est store abortionFullControl

est restore abortionFullControl
margins, at(treatmentGroup=(1(1)5) partisanship=(3) female=(0(1)1))

est restore abortionFullControl
margins, dydx(female) at(treatmentGroup=(5) partisanship=(3))

est restore abortionFullControl
margins, dydx(partisanship) at(treatmentGroup=(5) female=(1))





